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Abstract 



Recently, the B factory experiments BABAR and Belle as well as the CDF 
collaboration found evidence for mixing in the D meson system. The current 
status (beginning of summer 2008) of the experimental results of mixing is 
summarized. 



1 Introduction 

The most surprising result of last year's spring conferences on particle physics was 
the report of the evidence for — mixing by both the BABAR and Belle collab- 
orations [D [2]- At the end of 2007, the CDF collaboration found evidence for 
mixing [5] in a different environment but in the same decay channel as J3ABAR. 

The first mixing results were obtained in the neutral kaon system [4], 50 years 
ago. Mixing in the B^ system [5] was established in 1987 and measurements of the 
mixing parameters in the Bg system were published in 2006 by the CDF and DO 
collaborations [6]. 

Within the Standard Model (SM), mainly the first two generations contribute to 
the mixing and the mixing parameters and CP violation are expected to be very 
small. The observation of the mixing completes the picture of quark mixing, since 
the system, in contrast to the other three neutral systems, involves down-type 
quarks in the mixing loop. On the other hand, non-Standard Model processes could 
enhance either the mixing or the CP violation or both. Therefore, the measurement 
of large mixing parameters in the D*^ system or sizable CP violation would be a strong 
indication for New Physics (NP). 

In this paper, we present an overview of mixing. After an introduction to the 
charm mixing phenomenology and analysis techniques, results of the mixing para- 
meters and CP violation as related to mixing are summarized. They are obtained from 
hadronic two-body, multi-body final states and from quantum correlated decays 
of the experiments BABAR, Belle, Cleo and CDF. Mixing results from semileptonic 

decays can be found elsewhere [H [8] . 
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2 Mixing formalism and notation 




The neutral mesons are created as eigenstates of the strong interaction and can be 
distinguished by an internal quantum number, e.g. charm. Due to the weak interac- 
tion, an initially defined state of or l^^) will evolve with time into a mixture of 
D° and D^. The time evolution can be described by an effective weak Hamiltonian 
in the time dependent Schrodinger equation 

where M and T are mass and decay width matrices. The solutions of the Schrodinger 
equation are the mass eigenstates 

\D^,,) ^ p\D') T q\^) , 

with p'^ + — 1. They are linear combinations of the flavor eigenstates D° and !D°. 
The mass eigenstates, \Di) and |I?2), propagate independently in time with their own 
lifetime ri^2 and mass Mi^2 

|Di,2(t)) = e-^(^i-^-^^i-2/2)*|Z)i,2(t = 0)) . 

The ratios x = 2(Mi - M2)/{Ti + T2) and y = {Ti - r2)/(ri + T2) are related 
to the difference in lifetime and mass of the mass eigenstates. These variables are 
referred to as mixing parameters and are the observables to be measured. 

The probability / to find the state from an initial state after a time t is 



7(L>° ^ := |(L'°|L>°(i))|2 = ^[cosh(n) + cos(n)] 



-rt 
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and the one to find a l^^^) is 

/(D°^D°;t):= \(D'\D^{t))\^ = ^ 



2 
P 



[cosh(n) - cos(n)] 



After a certain time, the opposite flavor component appears. Mixing will occur if 
either the mass difference x or the lifetime difference y of the two states is non-zero. 
Depending on the size of x and y, an oscillating behavior can be observed, e.g. as in 

the case of the Bs system. 

If CP violation is neglected the state \Di) (I-D2)) is CP-even (CP-odd). However, 
CP violation in D mixing can be parametrised in terms of the quantities = ^ 

and (f)f = arg(^), where Aj = {flHolD^) ( Af = {f\HD\D^)) is the amplitude of 
a (-0°) to decay into a flnal state / and Tio is the Hamiltonian of the decay. A 
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value of 7^ 1 would indicate CP violation in mixing. A non-zero value of 0/ would 
indicate CP violation in the interference between mixing and decay. 

There are two contributions to the charm mixing processes, the short range box 
contributions and a long range part with on- and off-shell intermediate hadronic 
states. In contrast to the other neutral systems, the box diagrams of the sys- 
tem involve loops of down type quarks. Due to the GIM mechanism and the CKM 
matrix suppression, the lowest-order short distance calculation gives tiny results 
(a^box = O{10~^) and ?/box = 0(1O~^) [10]). Theoretical predictions for the long- 
range contributions are very difficult, as the c quark is either not heavy enough or 
too heavy to be treated by the different theoretical frameworks. Evaluations in the 
Operator Product Expansion (OPE) framework have shown that the mass and life- 
time differences are enhanced with increasing order in the OPE [llj compared to an 
OPE quark level analysis [12]. This yields O ~ 10~^ for both x and y. 

In the SM, CP violation in the charm sector is expected to be small and to be 
below the sensitivity of the experiments. Any measurement of CP violation in 
mixing would be a strong indication of NP. 

Calculations within the SM have large uncertainties and, therefore, the ability to 
detect NP contributions is limited. Nevertheless, the mixing measurements allow 
for a restriction of the parameter space of NP models. Golowich et al. (summary 
in [9|) have explored which NP models yield sizable values for x and y. They addressed 
21 NP models from various areas and found that in case of 17 models tightened 
restrictions on the model parameter space can be placed [T3] . 

3 Experimental techniques 

The most significant measurement of the mixing parameters was performed by the 
BABAR collaboration in the decay — > Ktt [1]. Therefore, basic ideas and techniques 
of mixing analyses are explained using this decay channel. 

3.1 Flavor tagging 

In order to perform a mixing measurement, the inital state has to be prepared and 
tagging of the fiavor at production time is required. The standard technique in charm 
physics is the use of D* D^n decays 0. Here, the charge of the pion TTtag determines 
the flavor of the at production time. 

The flavor at decay time is determined by the flnal state particle properties. Lets 
consider the decay of an initial which decays without mixing as K~tt^. 
In case of mixing, the converts to a D° and decays as K^tt^. Therefore, 

the charge of the K determines the flavor at decay time and thus if mixing occurred. 

^The charge conjugate modes are included throughout this paper. 
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Events of the decay mode ^ K tt^ are classified as riglit-sign (RS) and events of 
the decay mode ^ 71 are classified as wrong- sign (WS). 

Beside the Cabbibo-favored (CF) processes mentioned above, the D° also decays 
doubly Cabbibo-suppressed (DCS) as — > K'^ii'. As a consequence, the WS 
event sample contains not only events which have undergone mixing followed by a 
CF decay, but also DCS decays. The rate of the DCS events relative to the CF 
events is suppressed by a factor tan^(6'c7), where 9c is the Cabbibo angle. The rate of 
events with mixing is about another factor hundred smaller. Moreover, both decay 
amplitudes interfere. 

3.2 Extracting mixing parameters 

While the DCS events decay just exponentially, the mixed decays should have a more 
complex time structure due to the mixing process. The time evolution of the WS 
decay rate Tws can be approximated by 

TwsW oc e" 

\ " in / 

where Ru is the rate of the DCS events, y' and x' denote the mixing parameters 
and rt is the time in units of the decay time. CP conservation and small mixing 
parameters are assumed. There are three contributions to Tws, the DCS decays 
(I), the interference of the DCS and mixed decays (II) and the mixed decays (III). 
Each has a different time dependence. Measuring the time dependence of the WS 
decay rate allows to determine the mixing parameters x and y. Because of the strong 
phase difference between the CF decay amplitude and the DCS decay amplitude the 
mixing parameters are only defined up to a phase factor 6Kn- Therefore, the mixing 
parameters are measured as x' = x sin + y cos Sk-k and y' = —x sin Sk-k + 
y cosSktt- The phase 6k-k depends on the decay mode and in case of multibody 
decays it may vary over phase space. 

The mixing rate Ru defined as the time integral over the term III in Eq. [T]is inde- 
pendent of the strong phase (x^ + j/^ = x'^ + y''^) and can be measured in semileptonic 
D° decays. 

3.3 Event selection 

The Q value of the production process D* D^ntag is about 6 MeV which leads 
to a narrow peak in the difference of the D* and the mass (Am = m(D''7rtag) — 
m{Kn)). Selecting data in the peak region of the Am distribution suppresses back- 
grounds very effectively. The mass reconstructed from identified kaon and pion 
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tracks m(i^7r) is required to be in the expected mass window and the momentum 
in the center-of-mass system (CMS) has to be larger than 2.5 GeV in order to remove 
mesons from B decays. The kaon and pion tracks are refit to originate from the 
same vertex and form a which is fit together with the slow pion vr^ag to a common 
vertex. This provides an event-wise measurement of the proper time r and the 
error of the proper time a,-. The typical average value is 240 /xm with a resolution of 
100 yum. 

In the plane of Am and m{K7T), the BABAR analysis selects 1,129,000 RS and 
64,000 WS candidates from a data sample of 384 fb~^. 



4 Mixing measurements 
4.1 DO ^ Ktt 

The RS and WS event candidates in the plane of Am and m{Kn) contain differ- 
ent contributions of signal and backgrounds. Both are described by probability 
density functions (PDF). Their parameters are determined simultaneously for the 
RS and WS data sample in an unbinned maximum likelihood fit with four vari- 
ables m(/i'7r). Am, r, cr,- per event. Figure [1] shows the projections in m{KTi) and 




m{KTT)[GeVlc^\ m{KTT)[GeV / c^] 

Figure 1: m{KTT) (bottom part) and Am(i^7r) (top part) for wrong-sign candidates 
(right) and right-sign candidates (left). The fitted PDF's are overlaid. The colored 
regions represent the different background components [1]. 
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in Am(i^7r) of the RS and WS candidates. The black points denote the data and 
the sohd hnes show the fitted PDF's. The white surface below the curve is the ex- 
tracted signal. The dominant background component (blue) originates from properly 
reconstructed mesons combined with a random slow pion vrtag. 

From the fitted signal yields, the WS branching fraction -Rws in the decay D'^ —>■ 
Kit is extracted. The BABAR measurement [H] of i?ws = (0.353 ± 0.008 ± 0.004)% 
agrees well with the Belle result |I5] of i?ws = (0.377 ± 0.008 ± 0.005)%, where the 
uncertainties are statistical and systematic, respectively. 

The WS signal contains mainly DCS events with a small fraction of mixing events. 
In order to separate the mixing signal, decay time information is used as described 
in Sec. 13. 2[ The lifetime and the lifetime resolution function is determined from 
the RS event sample. The WS event sample is fit according to Eq. [1] The result is 
shown as solid curve in Figure [2l The mixing parameters are measured to be Rj) = 




Figure 2: a) Projections of the proper time distribution of the WS candidates and of 
the fit results. The results of the fit allowing (not allowing) mixing is overlaid as solid 
(dotted) curve. b) The points represent the difference between the data and the 
no- mixing fit. The solid curve represents the difference between fits with and without 
mixing pLj. 



(0.303±0.016±0.01)%, y' = (0.97±0.44±0.31)% and x'^ = (-0.022±0.03±0.021)%. 
The dotted curve shows the fit with the assumption of no-mixing. The histogram in 
the lower part of Figure [2] displays the difference of the data and the no-mixing fit 
(dots), while the curve is the difference between the mixing and the no- mixing fit 
model. The mixing model describes the deviations seen in the residuals. Thus, the 
deviations of the points from zero can be accounted for by the mixing. The 
significance of the mixing signal is evaluated based on the change of the negative log- 
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likelihood with respect to the maximum. Figure [3] (left part) shows the confidence- 
level contours calculated from the change in log-likelihood (— 2Aln£) in the two 
dimensions of x'^ and y' considering statistical uncertainties only. The likelihood 
maximum is indicated as black dot and is located in an unphysical region (x'^ < 0). 
The most likely physical allowed value (x'^ = and y' = 6.4-10"^) has a log-likelihood 
— 2A In C of 0.7 units. The value of — 2A In C for no-mixing is 23.9 units. Including the 
systematic uncertainties, this corresponds to a significance of 3.9 standard deviations 
and thus, constitutes evidence for mixing. 

To search for CP violation, Eq. [1] is applied to the WS and D ° samples 
separately, fitting for the parameters Rj), x'^ and y' for decays (-I-) and D ° decays 
(-). Both sets of mixing parameters (?/'+ = (0.98 ± 0.64 ± 0.45)%, x'+^ = (-0.024 ± 
0.043 ± 0.03)%) and {y'- = (0.96 ± 0.61 ± 0.43)%, x'-^ = (-0.020 ± 0.041 ± 0.029)%) 
are fully compatible with each other and differ by more than three standard deviations 
from the no-mixing hypothesis. The values i?D = \/ Rb-^i) — (0-303 ± 0.16 ± 0.10)% 
and Ab = {R^ - Rb)/{Rd + Rd) = (-2-1 ± 5.2 ± 1.5)% are calculated. Ab is fully 
compatible with zero. Both sets of mixing parameters do not differ. Therefore, no 
evidence for CP violation is observed. 
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Figure 3: (Left) Comparison of the mixing parameters x'^ and y' obtained by BABAR 
(black). Belle (red) and CDF II (green). The dark curves are the 1 — 3 standard de- 
viation likelihood contours from BABAR, the red curve is the two standard deviation 
contour from Belle (stat. uncertainties only) [Il[l5l[3]. (Right) Ratio of WS and RS 
decays as function of normalized proper time as obtained by the CDF collabora- 
tion. The dashed (dotted) curve is the fit with the mixing (no-mixing) hypothesis [3j. 



At the end of the year 2007, the CDF II collaboration published evidence for 
mixing in D*^ — > Kii decays using an integrated luminosity of 1.5 fb~^ in pp collisions 
[3j. The measurement exploits the time dependence of the number of WS and RS 

decays. In the range from 0.75 to 10 units of the decay time (Ft), (12.7 ± 
0.3) ■ 10^ WS and (3.044 ± 0.002) ■ 10^ RS signal events are selected. A least-squares 
parabolic fit of Eq. [T]to i? = A^ws/A^rs in 20 bins of Tt as shown in the right part 
of Figure [3] determines the mixing parameters i?D = (0.304 ± 0.055)%, y' = (0.85 ± 
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0.76)% and x'^ = (—0.012 ± 0.035)%. The no-mixing hypothesis is indicated as 
dotted hne. Despite the different production environment and analysis techniques, 
the agreement with the BABAR mixing resuhs is astonishing (see green dot in the left 
part of Figure [3]). Bayesian probability contours in the x'^ — y' plane are calculated. 
The data are inconsistent with the no-mixing hypothesis with a probability equivalent 
to 3.8 standard deviations. 

The Belle measurement of mixing in Kit published in 2006 is based on an 

integrated luminosity of 400 fb~^ [15j. Mixing parameters of i?D = (0.364 ± 0.017)%, 
y' = (0.06lg:^^)% and x'^ = (0.018t[}:i^)% are obtained. In Figure [3] the results 
are compared to the BABAR and CDF measurements. Belle excludes the no-mixing 
hypothesis with a significance of two standard deviations considering statistical un- 
certainties only. 

4.2 Strong phase 

The measurements of the mixing parameters x and y are only defined up to a strong 
phase Sktt (amplitude ratio of the CF to the DCS decays in —>■ Kn). The CLEO 
collaboration recently published a measurement of Sktt using quantum correlated 
pair^ which were produced in 281 pb ^ of e~^e collisions on the \Ef(3770) 
resonance [16]. Two general decay classes are considered, single tags and double 
tags. For the single tags, one D° is reconstructed independently of the other. This 
class provides uncorrelated decay information. In case of the double tags, both 
mesons are reconstructed, and they decay correlated. The final states can be hadronic 
or semileptonic. Depending on the final state, different enhancement factors to the 
ratio of the correlated and uncorrelated Z)° decay rates apply. The enhancement 
factors are functions of the mixing parameters and the strong phase. Therefore, the 
measurement of the time integrated yields of the correlated and uncorrelated 
decays allows to extract the mixing parameters and the strong phase by a fitting 
procedure. Using external branching fraction measurements in the fit, CLEO obtains 
cos{5ktt) = 1.03l!o;i7±0.06. Including additional external measurements of the mixing 
parameters in the fit provides an alternate measurement of cos{6Kn) = 1.10±0.35± 
0.06 and allows for a determination of x ■ sin (^Xtt) = (4.4l?] ± 0.29) ■ 10-3 and 
6k. = (22 tllt.lr. 

With these measurements, CLEO established a new technique of time-independent 
measurements of mixing parameters and the first measurement of the strong phase. 

4.3 ycp from lifetime measurements 

The decay time r of mesons (+) and mesons (-) decaying into final states of 
specific CP (such as K~K^ and tt^tt^) can be considered to first order as exponential 

2CLE0 uses the inverse amplitude ratio for the definition of the strong phase Sktt- 
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with a small correction term that depends on the mixing parameters: 



1 + COS0/ =F 2; sin0/) ' 

where is the lifetime of the CF decay Kir. The lifetimes can be combined 

into the quantities ycp and AY: 



T 







I/CP = — - 1 = ; , - 1 = ' - 1 AY 

with Ar = (r^ — r^)/(r^ + r^). In the limit of CP conservatioiifl, ycp = ±1/ and 
AY = 0. 

Belle found significantly different decay time distributions for decays to the CP- 
eigenstates K^K~ and tt+tt^ compared to the one to the CP-mixed state K~tt^ and 
measured ?/cp = (1.31 ± 0.32 ± 0.25)% with a significance of 3.2 standard deviations 
including systematic uncertainties f2]. The dataset corresponds to 540 fb~ . The 
measured lifetime asymmetry parameter A^ = —(0.01 ± 0.3 ± 0.15) • 10~^ reveals no 
evidence for CP violation. 

At the end of 2007, BABAR published a measurement of ycp from the lifetime 
of the three decay modes K^K~, 7r+7r~ and Kn [17] using 384 fb~^ of data. 
A value of ycp = (1.24 ± 0.39 ± 0.13)% is obtained, which is evidence for mixing 
at the three standard deviation level. No indication for CP violation was found, 
as indicated by the value Ay = (—0.26 ± 0.36 ± 0.08)%. Combining this result 
with a previous untagged BABAR measurement [18] yields the combined measurement 
ycp = (1.03 d= 0.33 ±0.19)%. 

The charm subgroup of the Heavy Flavour Averaging Group (HFAG) provides 
combined values for mixing parameters [22]. Besides BABAR and Belle, results from 
E791, FOCUS and CLEO contribute to the average values ycp = (1.132±0.266)% and 
Ar = —(0.123 ± 0.248)%. The precision is dominated by the B factory experiments. 
The ycp measurement clearly indicates mixing at a lifetime difference which is 
about 1 % and the no-mixing case is excluded at 4.5 standard deviations. There is no 
indication for CP violation from the averaged asymmetry measurements. 



4.4 Time-dependent Dalitz analysis 

Up to now, we considered only two body decays of the D'^. CLEO pioneered a 
method to measure x and y from a time-dependent Dalitz analysis of the resonant 
substructure in the decay K^tt^-k^ |_19j. This method allows to measure the 

sign of X. 

^Thc sign depends on the CP eigenvalue. 
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Belle extends the analysis of the self-conjugate process — > Kgn^n^ to a dataset 
of 540 fb~^ [20j. The Dalitz plot is shown in the upper left distribution of Figure |H 
It is described by a model containing contributions of 18 different quasi two body 
decays, which interfere. The time-dependent decay amplitudes of the and D ° are 
functions of the Dalitz variables {m'^{K^7r~^) , m'^{K^TT^)) and the mixing parameters 
X and y. Therefore, x and y can be extracted from an unbinned maximum likelihood 
fit to the Dalitz variables and the measured decay time. The data points are shown 
in Figure m and the curves represent the fit results. Assuming negligible CP violation. 




(GeVV) m^l (GeV^/c") 



Figure 4: Dalitz plot distribution and the projections in case of data (points with 
error bars) and the fit result (curve) as obtained by the Belle collaboration [20]. 

X = (0.80 ± 0.29lE]:[J?;j]i^)% and y = (0.33 ± 0.24l°:°^l°:°^)% are obtained, where the 
uncertainties are statistical, experimental systematic and decay-model systematic. 
The no-mixing case is excluded at 2.2 standard deviations. 

Allowing for CP violation, the fit of the additional parameters \q/p\ and indicates 
no evidence for CP violation in mixing or interference between mixed and unmixed 
amplitudes. Since the fit parameters are consistent for both the and D " sample, 
there is also no evidence for direct CP violation. 

BABAR finds evidence for mixing using a time-dependent amplitude analysis of 
the decay D° K~^tt~tt'^ in a data sample of 384 fb~^ [21]. The decay contains WS 
and RS events (see Sec. 13.11) . The signal and background yields are extracted from a 
binned extended maximum likelihood fit to the Am and m{KTni^) distributions. The 
time-dependent relative WS decay rate is a function of the Dalitz variables (m|^+^_ , 
m|,+^o) and the mixing parameters in the form x^^^o = x cos 5 k-k-kO + Z/sin^ft-^TrO 
and y^^^o = y cos Sk-k-kO — x sin Sk-k-kO , where Sk-k-kO is the strong phase between the 
DCS and the CF amplitude of p~K^ . The phase 5ktt-kO is different from 



10 



6ktt and has to be measured elsewhere. From a time-dependent fit to the WS data, 
x'^^^o = (2.6ll[J:^^ ± 0.39)% and y'^^^, = (-O.OGtgH ± 0.34) are derived with a 
correlation of —0.75. The significance is equivalent to 3.2 standard deviations. 

5 Combined results 

The HFAG determined world average values of the mixing parameters x and y in 
a global fit which takes into account all the relevant data from the various experi- 
ments [22] . Most of the more recent results have been presented in this talk. They are 
dominated by the B factory measurements with significant contributions from CLEO- 
c and Tevatron. The no-mixing case is excluded at about seven standard deviations 
{x = (0.91 ±0.26)% and y = (0.73 ±0.18)%). The mass difference differs from zero by 
3 standard deviations, the lifetime difference deviates from zero by about 4.1 standard 
deviations. The fit also determines Rb = (0.3342 ± 0.0083)%, 6k^ = (21.6l{^;^)° and 
SxttttO = {'^^■^-ll'lY- The measurement of CLEO of the strong hadronic phase Sk-k 
discussed in 14.21 was not considered because external measurements of i?D and Rm 
entered in the fitting procedure. 

Another set of parameters was determined allowing for CP violation. The values 
obtained for the CP sensitive parameters Ad = —2.2 ± 2.5, \q/p\ = 0.86lo;i5 and 
(p = (— 9.6lg;5)° indicate no evidence for CP violation within the current sensitivity 
of the experiments. 

The y measurement from the D° lifetime (see Sec. 14.31) yields y > . Therefore, 
the \Di) as CP-even state lives shorter than the CP-odd state \D2). The sign of x is 
measured in the Dalitz analysis (see Sec. 14. 4|) as a; > 0. Therefore, the CP-even state 
\Di) is heavier than the CP-odd state ID2) • 

Preliminary Monte Carlo studies of the LHCb collaboration indicate that the 
statistical uncertainties of the mixing parameters x, y and ycp measured in a time- 
dependent WS analysis of Kn decays and a lifetime analysis of decays to 
K~K~^ and K^n^ in 10 fb"^ of data may decrease by a factor five |23j . 

In summary, three experiments found evidence for mixing measuring lifetime 
and mass differences at the level of 1 %. The combined results of all experiments 
exclude the no-mixing case at seven standard deviations. The measurements are 
compatible with the SM expectations. It seems likely that the mixing is dominated 
by long-distance processes, which are difficult to calculate. Therefore, identifying NP 
contributions from mixing alone is not easily possible. Neither a single experiment 
nor the combination of all results provide any hint for CP violation in D° mixing. 
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